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Abstract 

A mathematical model is developed to describe the functional response of parasi- 
toids that imperfectly avoid superparasitism. This new model contains older models 
as special cases. It is satisfactorily fitted to the measured response of a laboratory 
parasitoid-host system. 



DEFINITIONS 

Parasitoids are insects requiring a host for the development of their larvae. The 
free-living adult females search for hosts and oviposit in or near them. When the eggs 
hatch, the larvae feed on the hosts, causing at first little damage, but gradually attacking 
more vital parts and eventually killing the host. 

Most parasitoids belong to two orders (Diptera and Hymenoptera) and they make 
up about 14% of the known insect species. Some of them are specific while others 
attack a variety of hosts species. Hosts belong to almost all insect orders, eggs, larvae, 
and pupae being far more frequently attacked than adults. Other species than insects 
can also be parasitized (spiders, snails). 

Parasitoids are called solitary if one larva only can develop in a given host, 
gregarious if several larvae can develop in a same host. Solitary species (e.g. most 
Ichneumonidae) normally lay a single egg at each oviposition. Gregarious species 
(e.g. most Braconidae and Chalcidoidea) usually lay a batch of eggs corresponding to 
the host’s capacity (up to several hundreds, even thousands) (Berland 1976). 
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Superparasitism is the situation in which more individuals of a parasitoid 
species occur in a host than can survive. This usually happens when a same host is 
parasitized more than once. When superparasitism occurs, competition takes place 
between the developing larvae and the supernumerary ones are eliminated by various 
means (Salt 1961). It is therefore a selective advantage to avoid superparasitism. That 
many parasitoids have evolved such an ability to some extent is a well established fact, 
for solitary as well as for gregarious species (e.g. Podoler & Mendel 1977). 

The number of hosts attacked per parasitoid per unit time changes with variation 
in host density. This relationship is called functional response (Solomon 1949). 
Holling (1959) distinguished three response types. The type II, which is the most 
common, is the only one considered here. It shows an initial increase from zero, with a 
gradually decreasing slope, tending to an upper horizontal asymptote (see Fig. 2). 

Variation in parasitoid density also affects the number of hosts attacked, but this 
relationship is not studied here. The parasitoid density is assumed to be constant. 



THE PROBLEM 

Existing models of the functional response of parasitoids deal with two cases: 

1) parasitoids that do not discriminate between healthy hosts and already parasitized 
hosts, i.e. parasitoids with no avoidance of superparasitism; 

2) parasitoids that recognize (and avoid) parasitized hosts both instantly and perfectly. 
In such a case, the hosts can be considered as being absent from the population and 
the parasitoids respond like predators, which remove the prey from the population. 
Models for these two cases were given by Rogers (1972). 

The model developed in this work deals with the more realistic case of parasitoids 
that spend some time inspecting every host encountered and that may fail to reject 
parasitized hosts. Such a model should contain models (1) and (2) as limit cases. 

The only previous attempt to build a model of this kind (Rogers & Hassell 1974) 
unfortunately suffers internal inconsistencies. 



THE MODEL 

Given P parasitoids and N hosts interacting for a time T, the model must evaluate 

the average number of hosts attacked per parasitoid per unit time, G, as a function 

of N and T. Let us make the following assumptions: 

1. The parasitoids search randomly, with a constant searching efficiency a (probability 
of a given parasitoid to find a given host). 

2. When finding a host, the parasitoids spend a time T e examining whether it is healthy 
or already parasitized. 

3. If the host is healthy, the parasitoid oviposits on it. 

4. If the host is parasitized, the parasitoid has a probability p of making a “mistake” 
and ovipositing on it. 

5. During and after each oviposition the parasitoid does not search for a time T 0 
(time for ovipositing per se, time for resting, etc.). 
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These assumptions are summarized on a diagram (Fig. 1). They yield the following 
expression for the functional response G (the mathematical proof is given in the 
Appendix) : 




where T h = T e + T 0 is the time spent “handling” a healthy host and T p = T e + pT 0 
is the average time spent “handling” a parasitized host. 

Equation (M) is an implicit equation that must be solved by numerical iteration. 
Interestingly, the number of parameters needed is one less than the number that were 
introduced in the hypotheses. The number of hosts attacked per parasitoid per unit 
time is completely determined by 7, A r , P, a, T h , and T p . The last three parameters 
characterize the parasitoid species (in connexion with a given host species). 




Hg. 1. 

Conceptual diagram of the model for the interaction of P parasitoids and N hosts for a time T. 



AN EXAMPLE 

The Chalcidoid gregarious parasitoid Nasonia vitripennis (Walker) attacks the pupae 
of the house fly Musca domestica L. The functional response of this system in laboratory 
conditions was measured by Edwards (1961). 

The new model (M) was fitted to these data, using a standard non-linear regression 
method (Fig. 2). With P — 8 parasitoids and T — 24 hours, the following values were 
obtained for the parameters: a = 0.44 day -1 , T h = 0.20 day, and T p = 0.0S day. 
According to these estimates the parasitoid spends therefore 2.5 times more time 
handling a healthy host than a parasitized host. 




HOSTS ATT ACKED/DAY/PARASITOID 
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Fitting alternative models to the same data yielded the following estimates: 

1) “Random parasite” model (assumption T p = T h , see Appendix, eq. (12)): a = 1.03 
day' 1 and T h = 0.19 day. 

2) “Random predator” model (assumption T p = 0, see Appendix, eq. (13)): a = 0.32 
day' 1 and T h = 0.21 day. 




HOSTS/CONTAINER 



Fig. 2. 

The functional response of Nasonia vitripennis (Hymenoptera : Pteromalidae) 
attacking house fly pupae, with eight parasitoids searching simultaneously (Edwards 1961). 
The curve is the least squares fit of model (M). 



The quality of the fit to the data is equally good with all three models but the 
estimated values of the searching efficiency a are rather different. As could be expected, 
the estimate of the new model is intermediate between those of the other two extreme 
models. It is not possible to select one model against the others with purely statistical 
arguments, but biological knowledge of the parasitoid’s behaviour is necessary. 

Such information is given by Wylie (1965): Nasonia females always drill first 
through Muscat puparium, and after piercing the enclosed pupa with the ovipositor 
usually (but not always) withdraw rapidly from an already parasitized pupa. Thus, 
there is a definite discriminatory capacity, which is neither immediate nor perfect. The 
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new model should then be preferred in this case. But to reliably validate this model, 
estimated values of the parameters (especially the searching efficiency a) should be 
compared to direct observational measurements, which are not available. 



A mathematical model is presented, which deals with the problem of insect para- 
sitoids that imperfectly discriminate between healthy and parasitized hosts. The aspect 
examined is the functional relationship between host density and the number of hosts 
attacked per parasitoid per unit time. 

On the basis of a few behavioural assumptions, the model is at first built concep- 
tually and then translated into mathematical equations. The model fully characterizes 
a given parasitoid with three parameters: searching efficiency, time spent handling a 
healthy host, and time spent handling a parasitized host. 

The model is fitted to the reported functional response of an experimental para- 
sitoid-host system, and the parameters are identified. The comparison of these estimates 
to those obtained from other existing models is then briefly discussed. Although con- 
ceptually more satisfactory, the new model needs additional experimental data to be 
reliably validated. 



Summary 



Appendix 



The derivation of expression (M) is given here. 

The number of hosts attacked during the time period T must be calculated. 
N h (t) and N p (t) vary with time t but 



N h (t) + N p (t) = N = const. 



(i) 



Let us first describe the activity of one parasitoid. 
Let dt be a short time interval. Let 



time spent 
searching 




( time spent 
examining 
hosts 




time spent 
ovipositing 



dt = d t s + d t e + d t 0 



( 2 ) 



Number of encounters with healthy hosts during dt 

d E h = aN h (t) d t s 



( 3 ) 



Number of encounters with parasitized hosts during dt : 

d E p = aN p (t)dt s 

We have, by definition of T e , T 0 and p: 



( 4 ) 



d t e = T e (dE h +d£ p ) 
dt 0 = T 0 (dE h +pdEp) 



( 5 ) 

( 6 ) 



892 



R. ARDITI 



Inserting (3)-(6) into (2) yields 

df = dt s { 1 + a (T e + T 0 ) N h (t) + a (T e +pT 0 )N p (f ) } (7) 



Let us define 




+ T„ 


( 8 ) 


+ PT 0 


( 9 ) 



The change in the number of healthy hosts during d t, due to attacks by all para- 
sitoids, is 



d N h = —PdE h = —aPN h (t) dt s 



With (7)-(9) and (1) we obtain the following equation for the dynamics of the 
population of healthy hosts : 

d Nh = -aPNj, (0 (10) 

d t 1 +aT p N +(T h -T p )N h (t) 

This first order nonlinear differential equation must be integrated between t = 0 
and t — T. Assume that all hosts are healthy at t = 0, i.e. N h (0) = N. Let = N p ( T ) 
be the total number of hosts attacked. 

The integration of eq. (10) yields the following transcendental equation for N 

N "- N ( i -' xp (^ (r ' <T ‘- r ' )1 ? ) )) (n) 

By definition, the functional response is given by G = Npl(PT), where G is the 
number of hosts attacked per parasitoid per unit time. Writing eq. (11) for G instead 
of Np yields expression (M). 

The existing models of Rogers (1972) appear as special cases of eq. (11): 

1) T p = T h , i.e. p — 1 (parasitoids do not recognize parasitized hosts). 

- N (‘ - exp (12) 



This is the “random parasite equation”. 



2) T p = 0, i.e. T e = 0 and p = 0 (parasitoids recognize parasitized hosts immediately 
and without failure). 



N T P = N 




a (TP-T h N 




(13) 



This is the “random predator equation”. 
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